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This paper develops the properties of the Erlang loss function, B(N, a), 
used in telephone traffic engineering. The extension to a transcendental 
function of two complex variables is constructed, thus permitting the 
methods of complex analysis to be employed for the further study of its 
properties. Exact representations, Rodrigues formulas, and addition 
theorems are given both for the loss function and for the related Poisson- 
Charlier polynomials. Asymptotic formulas and approximations are 
developed for the loss function and also for its derivatives. A table of 
coefficients is included which, together with one of the asymptotic formulas, 
permits computation of B(N, a) by simple means even when the number 
of trunks, N, is very large. This same table is used to obtain dB(x, a)/dx. 

I. INTRODUCTION 

The Erlang loss function 

B w>*~m/kt\ (1) 

is fundamental to the study of telephone trunking problems. A. K. 
Erlang 1 used B{N, a) to express the probability that a call, which is a 
member of a Poisson stream of parameter a, arriving at a group of N 
telephone trunks will be rejected. Later studies of trunking problems 
have shown the desirability of enlarging the scope of applications of 
the loss function. For example, the consideration of trunk groups with 
nonintegral number of trunks arises in determining the equivalent 
number of trunks in Wilkinson's "equivalent random method." 2 
Methods for accomplishing the computation by interpolation are given 
by Rapp 3 while continued fraction procedures for accurate computa- 
tion are given by Levy-Soussan 4 and Burke. 5 Derivatives with respect 
to N and a arise in optimal trunk group size apportionment problems. 
See, for example, Akimaru and Nishimura 6,7 who studied such models 
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and prepared tables of derivatives. In some investigations, rapid and 
accurate approximate computations of B(N, a) for very large trunk 
groups are needed. This occurred in the study of certain satellite 
telephonic communication problems. 89 The need thus arises of en- 
larging the definition of B(N, a) as given in (1). Of course, that is done 
implicitly in the above investigations. It has been customary to extend 
the definition of B(N, a) by use of an integral formula (Theorem 3) 
ascribed to Fortet. This integral formula is used in (23) to define a 
transcendent, B (z, a) , for complex z and a. The extension to the complex 
plane in both z and a permits the powerful methods of complex analysis 
to be applied for obtaining exact, asymptotic, and approximate 
representations. 

It is the purpose of this paper to provide an investigation into the 
properties of B{z, a) with the object of generalizing known results, 
obtaining new results, and presenting practical methods for application 
to the class of problems outlined above. 

Part II derives exact relations satisfied by B(z, a). Similar relations 
for the related Poisson-Charlier polynomials, G } (z, a), are derived in 
the appendix. These relations provide efficient means for exact com- 
putation; thus, Theorems 1 and 2 constitute a practical method of 
computing B(N, a) to a prescribed accuracy for isolated computa- 
tions. Similarly, the use of Theorem 5 enables one to compute B(z, a) 
even for nonintegral number of trunks. Theorem 6 may be similarly 
employed. The relationship of B(z, a) and Gj(z, a) to Whittaker func- 
tions as given in Theorems 7 and 24 is the key for linking up these 
functions with the more well-known functions of applied mathematics, 
i.e., hypergeometric functions and Laguerre polynomials. The Rod- 
rigues Theorems 8 and 22 are useful for the evaluation of integrals 
of the form 

fp r >(a)a- 1 e-"B(N, a)~Ha, // (r) (a)o««— <?,(*, a)da (2) 

and, as in the case of Theorem 22, for obtaining an integral representa- 
tion. The addition Theorems 9, 10, 26, and 27 yield group-theoretic 
structure information which is useful for simplifying formulas contain- 
ing these functions, and for the evaluation of integrals. The evaluation 
of an integral, by means of generating functions and Theorem 10, was 
done in Part IV to obtain ultimately an approximate formula for 
dB(x, a)/dx. A general use of the exact relations is to serve as a spring- 
board for asymptotic and approximate results and also for their error 
estimations. This is well illustrated in Part III of the paper. 
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The asymptotic expansions of Part III are also representations of 
B(z, a) but, unlike those of Part II, when used as approximate formulas 
for computation they cannot yield results of arbitrarily high accuracy, 
i.e., the accuracy depends on specific values of parameters. Theorem 
11 is particularly useful for computation when \z/a\ is small. It may 
be used for the computation of B (z,a) for fractional number of trunks 
by computing B(z, a) for < z < 1 and then using the recurrence 
formula of Theorem 4. Theorem 11 includes well-known asymptotic 
results, e.g., 



(- |, oV = vW(l - erf Va) ~ 1 - \ o~» + ^ cr* 

a~ 3 -\ , a -* * , (3) 



B 
\ v ' / 

_ L3J) 

2 3 

B(-l, a)" 1 = - ae a Ei{-a) ~ 1 - a" 1 + 2 ! er 2 

- 3 ! a- 3 H , a -> oo . 

An undesirable feature of many methods of computing B(x, a) is the 
dependence of the computational effort, e.g., time of computation, on 
the value of x ; thus, the larger the value of x the greater the computa- 
tional effort. Theorem 14 overcomes this defect; the computational 
effort is independent of the size of x. Theorem 14 is easily usable even 
with a desk machine regardless of how large x is. The accuracy, how- 
ever, depends on x and a parameter c. For fixed c the accuracy improves 
with increasing x. When x is fixed, the accuracy deteriorates when c 
is large and negative but greatly improves as c is increased. To facilitate 
the use of Theorem 14, Table I gives required coefficients, namely, 
cio(c), cii(c), ai{c). To use the table, one computes 

c = — r ~ , (4) 

\x 

then 

B(x, a)" 1 ^ a (c)Vx + oi(c) 4- ^. (5) 

Possibly, one should comment that the range of values of x, c for 
which (5) is accurate is not as important as the fact that it is accurate 
over a wide range of values of B(x, a), that is, values encompassing 
the ranges of most applications. For quantitative limitations, see Fig. 
1. A method of obtaining dB(x, a)/dx based on Theorem 14 is given 
in Part IV. This uses the formula 
dB (x, a) ^ B(x, a) 2 / _oA 
9x 2vx \ -c / 



.• _l_ r, I -r 1 

(6) 



a B(x, a) (-- 1 + B{x,a) 



2x 
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Table I — Coefficients for evaluation of B(x, a) and 3B(x, a)/dx 



c* 


a 


ai 


a 2 


c* 


a 


ai 


a 2 


-3.0 


225.3 


2032 


13726 


0.6 


0.8230 


0.7274 


0.1011 


-2.9 


167.7 


1367 


8536 


0.7 


0.7749 


0.7414 


0.0985 


-2.8 


126.0 


925.4 


5334 


0.8 


0.7313 


0.7552 


0.0954 


-2.7 


95.63 


630.5 


3348 


0.9 


0.6917 


0.7686 


0.0920 


-2.6 


73.28 


432.2 


2111 


1.0 


0.6557 


0.7814 


0.0883 


-2.5 


56.70 


298.0 


1336 


1.1 


0.6227 


0.7937 


0.0845 


-2.4 


44.29 


206.7 


848.1 


1.2 


0.5926 


0.8053 


0.0806 


-2.3 


34.92 


144.1 


540.2 


1.3 


0.5649 


0.8163 


0.0767 


-2.2 


27.80 


100.9 


345.0 


1.4 


0.5394 


0.8267 


0.0729 


-2.1 


22.33 


71.07 


220.7 


1.5 


0.5158 


0.8364 


0.0691 


-2.0 


18.10 


50.27 


141.4 


1.6 


0.4940 


0.8455 


0.0654 


-1.9 


14.80 


35.71 


90.70 


1.7 


0.4739 


0.8540 


0.0619 


-1.8 


12.21 


25.49 


58.17 


1.8 


0.4551 


0.8619 


0.0585 


-1.7 


10.16 


18.27 


37.28 


1.9 


0.4376 


0.8694 


0.0552 


-1.6 


8.521 


13.15 


23.86 


2.0 


0.4214 


0.8763 


0.0521 


-1.5 


7.205 


9.522 


15.23 


2.1 


0.4062 


0.8828 


0.0492 


-1.4 


6.139 


6.936 


9.692 


2.2 


0.3919 


0.8889 


0.0464 


-1.3 


5.271 


5.090 


6.141 


2.3 


0.3786 


0.8946 


0.0438 


-1.2 


4.557 


3.772 


3.872 


2.4 


0.3661 


0.8999 


0.0413 


-1.1 


3.968 


2.830 


2.430 


2.5 


0.3543 


0.9049 


0.0390 


-1.0 


3.477 


2.159 


1.519 


2.6 


0.3432 


0.9095 


0.0368 


-0.9 


3.066 


1.682 


0.9486 


2.7 


0.3327 


0.9139 


0.0347 


-0.8 


2.721 


1.344 


0.5960 


2.8 


0.3228 


0.9179 


0.0328 


-0.7 


2.428 


1.108 


0.3816 


2.9 


0.3134 


0.9218 


0.0309 


-0.6 


2.178 


0.9435 


0.2540 


3.0 


0.3046 


0.9254 


0.0292 


-0.5 


1.964 


0.8318 


0.1804 


3.1 


0.2962 


0.9287 


0.0276 


-0.4 


1.780 


0.7580 


0.1398 


3.2 


0.2882 


0.9319 


0.0261 


-0.3 


1.620 


0.7112 


0.1187 


3.3 


0.2806 


0.9349 


0.0247 


-0.2 


1.481 


0.6840 


0.1089 


3.4 


0.2734 


0.9377 


0.0234 


-0.1 


1.360 


0.6705 


0.1052 


3.5 


0.2666 


0.9403 


0.0222 





1.253 


0.6667 


0.1044 


3.6 


0.2600 


0.9428 


0.0210 


0.1 


1.159 


0.6696 


0.1048 


3.7 


0.2538 


0.9451 


0.0199 


0.2 


1.076 


0.6771 


0.1052 


3.8 


0.2478 


0.9473 


0.0189 


0.3 


1.002 


0.6877 


0.1052 


3.9 


0.2421 


0.9494 


0.0179 


0.4 


0.9357 


0.7000 


0.1045 


4.0 


0.2367 


0.9514 


0.0170 


0.5 


0.8764 


0.7135 


0.1031 











n rr 

Calculate c = — — , then 



B(x, a) -1 ^ a Vx + ai + a 2 /Vx, 
dB(x, a) B(x,aY( a 2 \ 



dx 



2Vx 



^ *<*■<"{; -!+*<*,<•>}• 



Standardized offered load. 



It is appropriate to mention, at this point, another method of approxi- 
mating B(x, a) by means of a formula whose computational effort 
is also independent of x and which, similarly, is applicable over a wide 
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POINTS IN UNSHADED AREA 
YIELD ERROR < 5% 




-2.0 -1.5 -1.0 

C = STANDARDIZED OFFERED LOAD 

Fig. 1 — Five-percent-error contour. 



-0.5 



range of values of B(x, a). This method is described in Refs. 10 and 11. 
A comparison of this method with that of Theorems 1 and 2 is given 
at length in a report by S. Miller. 9 

Derivatives and inequalities on derivatives arc given in Part IV. 
Theorem 15 extends the well-known derivative formula for B(x, a) 
with respect to the real variable a. Theorem 17 provides an accurate 
approximation for dB(x, a)/ dr. Empirically, the accuracy seems to hold 
to four significant figures or better over a very wide range of values 
of x and a. Of significance is the corollary which shows that the ap- 
proximate value obtained is always too small. If a quick appraisal of 
the derivative is desired, Theorem 18 may be used. The logarithmic 
convexity properties of B(x, a) given in Theorem 19 provide the useful 
bounds of the corollary on the second derivatives. Also an application 
is given to the logarithmic interpolation of Theorem 20. This is very 
useful when, for example, one wishes to compute B(x, a) for x be- 
tween consecutive integers, say N, N -f- 1, and for which B(N, a), 
B(N + 1, a) are known. An extension of this idea is provided by 
Theorem 21 which permits accurate computation of B(x, a). 

It may be remarked that generally relations, representations, and 
asymptotics for B(z, a)~ l are simpler in structure than those for 
B{z, a) and may provide greater numerical accuracy in computations. 
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II. CONVERGENT REPRESENTATIONS 

The study of telephone trunking problems, whether in equilibrium or 
transient condition, or even nonstationary, 12 engenders the Erlang 
loss function, B(N, a), which initially arises in the form 13 

B(N, a) = j£j /e j y N ^ (integral), a > 0. (7) 

For these reasons and for the purposes of studying certain forms arising 
in queuing theory related to B(N, a) and also for the facilitation of 
numerical evaluation, it is useful to represent the loss function in 
diverse ways. 

The numerical computation of B(N, a) as given in (7) is awkward 
when a and N are large since then both numerator and denominator 
are large. A form well adapted to numerical work is 

& (8) 

# a.) = i f tfU) = jv"(JV - 1) ■ • • (N - j + 1) (j > 0), 

which follows from 

B(N, a)"' =| o f/a- . 1^,.- - £*«„-, (9) 

A modified form of (8) is given in Theorem 1. 
Theorem 1 : 

B(N, a)- 1 - E N^ar* + N^ar'B(N - v, a)~\ v ^ 0. 



E 

Proof: Since 



iVO+o = N lr) (N - „)«>, (10) 

one has, from 

B(JV, a)- 1 = E tf«>«r' + E AT^a->, (11) 

E tf M) er' - E N u+r) a- j -' = N M ar' E" (# - «0 (i) a _i 

/-> i-o J-° 

= N (v) a~'B{N - v, a)" 1 . (12) 

The formula of the theorem follows from (11) and (12). 
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Corollary: The case v = 1 implies the known 1 * difference equation 

B(N,a) = ^ • 

1 + 



aB(N - 1, a) 
R. Franks suggested using the value of B V (N, a) denned by 

B„(N, a) = 1 /"L N<»ari (13) 

/ j-0 

to approximate B(N, a) in which, for any small number r? > 0, the 
index v is chosen so that 

N M a- ^ ,,. (14) 

Theorem 2 bounds the error of the method. 

Theorem 2 : 

3,(N, o)(l - v) £ B(N, a) ^ B V (N, a). 

Proof: From Theorem 1 one has 

-= = B(N, a) ^ B V (N, a). (15) 

B V (N, a)" 1 + N^a"B{N - v, a)- 1 

Thus 

B(N, a) m 1 (16) 

B y (N, a) 1 + N^a->B,(N, a)B(N - v, a)" 1 ' 

Since N M aT' is strictly monotone increasing as a function of N, (8) 
shows that 

B{N + 1, a) < B(N, a) (17) 

for all N ^ ; thus 

Ml a > > - l - , (18) 

B,(N, a) ~ 1 + vB.(N, a)B(N, a)~ l 

and hence 

g (Ar ' a) > 1 - ,. (19) 

B„(N, a) ~ 

The theorem follows from (15) and (19). 
An integral representation, ascribed to Fortet, 15 may be obtained for 

B(N, a). 



Theorem 3: 

B 



(N^)-* = aj\-°»(l + y) N dy. 
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Proof: From the Eulerian integral 

J" e-"y l dy = T(l + l)a-'-\ I > - 1, (20) 

one obtains 

N l »ar j = a ( ■) J* e- av y j dy. (21) 

Use of (8) now yields 

B(N, a)" 1 = aT er*» £ ( N ) y j dy = af e-°»(l + y) N dy. (22) 
Jo y=o \ J / Jo 

The integral representation now permits extending B{N, a) into 
the complex plane with respect to both N and a. One defines 

B(z, a)- 1 = aT e-*»(l + yYdy (23) 

Jo 

in which 2, a may both be complex. Clearly, B(z, a)~ l is an entire func- 
tion of z for Re a > (Re designates "real part"). The symbols N, a 
will be used for nonnegative integers and positive reals, respectively. 
A generalization of Theorem 1 is given in Theorem 4. 

Theorem 4 '■ 

B(z, a)' 1 = E z^a~ j + Z<>>cT'B(z - v, a)~\ Re a > 0. 



;=() 



Proof: Integration by parts of (23). 

It is of interest to investigate the relationship of B (z, a) to the func- 
tion 

+ (z ' a) = e ~ a WTTy (24) 

which is an extension of the Poisson distribution function, \f/{N, a), 
with parameter a. The function if/(N, a) is a good approximation to 
B{N, a) when o is much less than N. Exact relations between B (z, a) 
and ^(z, a) are given in Theorems 5 and 6. These relations provide 
convenient means of calculation of B(z, a) for general z, a; e.g., in 
trunk group blocking problems when a nonintegral number of trunks 
is considered. 

Theorem 5 : 

B (v«)-' = * (Z - a) " " %i (> + \y~ (z + .r 

532 THE BELL SYSTEM TECHNICAL JOURNAL, MARCH 1974 



The series converges uniformly everywhere in Re z > — 1, Re a > 0. 
Proof: Let u = 1 + y in (23), then 

B(z, a) -1 = ae a J e- au u z du; (25) 

hence, 

B{z, a)~ l = ae a f e~ au u z du - a I e a(l - u) u z du, (26) 

Jo Jo 

and 

B(z, a)' 1 = + (z, a)" 1 - a f e a(l ~^u z du. (27) 

To exhibit the integral in (27) as an inverse factorial series, consider 
the beta function integral 

which yields the special case (s ^ integral) 

/„' M ' (1 " u) ' du = (z + l)..°l + l + s) ■ (29) 

Use of the expansion 

e a(l-«) = £ £L" (i _ M )« (30) 

s =0 S ! 

in (27) and subsequent use of (29) yield the result of the theorem. 
The Mittag-Lefner expansion for the integral of (27) leads to 

Theorem 6: 

Conditions of convergence are the same as in Theorem 5. 
Proof: The expansion 

e— = f (-1)*^ (3D 

s=0 s ! 

used in (27) leads immediately to the required result. 

Whittaker functions, 16 W k , m (z), play a useful role in the discussion 
of B (z, a) and of Poisson-Charlier polynomials to be introduced later. 
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They may be introduced by 

*...» = m : ; + m) / o - (i + -: ) f-** (3^ 

Re (| - fc + m) > 0, |arg z| < 7T. 

Theorem 7 : 

J5(z, a)" 1 = e- z/2 e ia ^ z /2,(,+i)/ 2 (a). 

Proo/: Let t = ay in (23), then 

B(z, a)" 1 = f" e-' (l + ^ )**• ( 33 ) 

The required result follows on comparison with (32). 

A Rodrigues type of relation for B(N, a) -1 may be obtained from 
Theorem 3. 



B(N + M, a)" 1 = {-\) M ae" -^[e-o^B (iV, a)" 1 ]. 



Theorem 8 : 

Proo/: From Theorem 3, one has 

e-°a-'B(N, a)~ l = T «-»+«'> (1 + y) N dy; (34) 

hence, 

(_l)jr ^ [g-a-iBQf, a)- 1 ] = /"" *-*+»> (1 + y)*+*d». (35) 
aa M Jo 

The formula for B(iV + M, a) -1 now follows on multiplication by ae a . 
Corollary : 

tff\T „\— 1 — ( \\N no a 

da 1 



B(N, a)- 1 - (-D^ae-^rCe-a- 1 ]. 



Additional formulas for B(z, a)" 1 (as a function of a) provide con- 
venient means of computation for values of a near some fixed point. 
Two such formulas are given in Theorems 9 and 10. 



Theorem 9 : 
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Proof: The function S„{a) given by 

S n (a) = ± - x (36) 

r =0 T ! 

is an Appell polynomial, that is, 

^ = S_,(«). (37) 

Thus the Taylor expansion for S n (a + t) can be written in the form 

S n (a + t) = t -.S n -,(<*)t>. (38) 

One obtains from (7) 

5(n, a )-' = n! a -»5nW (39) 

and hence the theorem follows on substitution into (38). 
Theorem 10: 

B(z, a + 0* 1 = (l+-)e<it t^ B(z + v, a)~H', 

\ a) „=o v\ 

Re a > 0, 1 1 1 < Re a. 

Proof: Let 

l(z, a) = e iz '- a a- l B(z, a)' 1 , (40) 

then, from (23), 

d d f* 

-j-l{z,a) = -7- e"' / e- a < 1+ «)(l + y)'dy, 

= - e izr [ e- a < 1+ «')(l + y) t+1 dy, 

Jo 



= K* + 1, a), 


(41) 


and hence 




d" 


(42) 


Thus, by Taylor's formula, 




/(2,a + = ± -.l(z+v,a). 

v =0 v ! 


(43) 



Substitution of (40) into (43) yields the required result. One has 
l(z 4- v, a) = ««■+'>« /"" e— «+»>(l + y)* + '(fy, (44) 
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hence the terms of (43) are 0[(*/Keo) '»■■•]• Tne stated convergence 
criterion now follows. 

III. ASYMPTOTIC EXPANSIONS 

Particularly simple and convenient forms for theoretical and numer- 
ical applications may be obtained by examining asymptotic expansions. 

Theorem 1 1 : 

00 

BM)" 1 ^!^"'! a— > ■», |arga| <t. 

Proof: The asymptotic expansion for W k , m (z) is 16 
W k , m (z) ~ e-*'z k 



X 



[m 2 - (k - £) 2 ][m 2 - (fc- j) 2 ]- • -[m 2 - (fc- y+*)»] 



1+ ^" v!z 



•=i 



z — * oo , I arg 2 1 < it. (45) 



Substitution of the parameter values given by Theorem 7 establishes 
the result. 

It should be remarked that the error, when using the partial sum 
Y. l t=l z< - v) or' to approximate Biz, a) -1 , does not exceed |a| \z {k) a- k \/ 
Re a provided Re z ^ k, Re a > 0. This follows directly from Theorem 
4 and (23). 

For large z, one has 

Theorem 12: 

B(z, «)-> - Hz, a)" 1 " t (0+ !).% + ,) . 

z — * oo , | arg z | < x/2, uniformly in any bounded region of the a-plane 
for which Re a > 0. 

Proof: The representation of Theorem 5 is used. One must show 



B(z, «)-i - Hz, «)-i - L 



= o ( ( z + i)- a " (* + *) )> (46) 

that is, 

lim £ 7—; ,"?" p , v = 0. (47) 

«-. a>„ (z + n + 1) • • • (z + s) 
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Let Re z = x. then one has 



1 



Let t; = x 4- n and / = s — n, then the dexter of (48) is 



Use of (29) and (30) on (49) yields 



• (48) 



(49) 



■r»(* + n + l)---(* + «) 



thus 



.r«(* + n +!)■•■ ^ + «) 



£|a| /" l e' a i (1 - u Vdw; (50) 

0, »->«>. (51) 



a 6 



» + l 



The theorem is proved. 

Useful asymptotic formulas are obtained when both a and z have 
infinite limits but approach infinity in a fixed ratio, that is, a = cz, c 
fixed. The cases c > 1, c = 1 are discussed by A. Descloux 17 for large 
real z. Theorem 13 generalizes the result for c > 1 to complex z and 
provides the structure of the coefficients for the complete expansion. 
The case c = 1 is obtained as a corollary to Theorem 14 where the 
result is also generalized to complex z. 



Theorem 13: 



B(z, cz)- 1 ~ £ q*r\ 

J=0 



z-+ oo, |argz| < 2 , c > 1, 

/ c d V c 

(Jl ~ \ c - 1 dc ) c - 1 " 

Proof:' One has, from (23) 

B(z, cz)~ l = cz I e-""(l + y)'dy. 



(52) 



* The author wishes to thank C. L. Mallows for this proof, which replaces a much 
longer proof originally supplied by the author. 
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Defining the function h{y) by 

h(y) = cy - In (1 + y), (53) 

one may write, since MO) =0, Ji(°°) = co , and Mi/) is monotonic 
increasing, 

^,«)--»/;«-"" . ( i ( +t)-i ,ft - (54) 

The factor c(l + y)/[c(l + y) - 1] is now expanded in powers of h 
as follows : 

,M+tA °S7i*- (55) 

c(l + y) — 1 j=o t! 

A theorem on Abelian asymptotics for Laplace transforms 18 and (54), 
(55) yield the asymptotic behavior of B (z, cz)~ l for «-*<», | arg z \ < ir/2 ; 
thus, 

B{z,cz)- l ~ iLgtr 1 - (56) 

1=0 

The coefficients gi may be evaluated as follows. Let 

w = c — In c, (57) 



and 
then 



fc( w ) = — £_ , (58) 

c — 1 



fe(u; + A) = fc[c(l + ?/) - In c(l + i/)] = c(1 C +y) _ t 



(59) 



Thus, Taylor expansion yields 



c 
One has 



(61) 



dw c — 1 dc' 
hence 



The following formula is obtained directly from Theorem 13. 



(62) 



n , c c 1 . 2c 2 + c 1 fR „s 

B % «)- ~ — , - j—^ - + ^^p 72 • (63) 

The evaluation and behavior of B(z, a) for a in a neighborhood of z 
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is often of interest; accordingly, the function B(z, z + cyfz) ' will be 
considered for z — > oo ; c is a fixed real number. 

Theorem 14'- There exists a representation of the form 
B(z, z + cyfz) ~ E aj {c)z-^-^\ 

;=0 



I ar g z I < o » c rea ^' 



in which 

ao(c) = e icl / e~ iu, c/M, 

2 1 2 _ 1 
S^S 3 



2 1 1 

ai(c) = - + -c 2 - ^ c 3 a (c), 



Proof: From (23), one has 

B(z, z + cy[z)~' = (z + cVz) r e -(^^)"(l + w)«dw, (64) 

Jo 



| arg z | < ^ 



2 

Let u = v/yfz, then 

B( 2 , 2 + cyfz)~ l = f M e-^'+^Hv, z)dv, 

/ v (65) 

h(v, z) = e*< 2 -^° ( 1 + ^ j (yTz + c). 

Let if be a positive constant, then, for | v | ^ K, h(v, z) clearly possesses 
an asymptotic development in yfz uniformly in v ; thus, 

h(v, z) ~ t bi(v, c)z-^>\ z^cc, (66) 

in which the coefficients b/(fl, c) are polynomials in u. In particular, 
6 (y, c) = 1, 
bi(», c) = \v 3 + c, (67) 



3 

1 

3 "' 1 ' 18 
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bt(v, c) = kCv 3 - 7» 4 + ia»'. 



Since 

e-(*»n-«>f,* e L(0, °°) (68) 

for each k > and any c, termwise integration of (66) leads to the 
required asymptotic expansion. Thus, letting 

aj(c) = f X *-«"■+">&*(», c)efo, (69) 

Jo 

one has 

B(«, 2 + cV^)- 1 ~ £ o;(c)2-"- 1 " 2 . (70) 

The formulas for a (c), ai(c), a 2 (c) stated in the theorem are obtained 
by evaluation of (69) using bj(v, c) as given in (67). 

Corollary : 

B(z ,w_ /« + 2 ■ A /z 

BV, V V 2 3 ^ 12 \2z' 

z-> co, |argz| < 2* 

Proof: The result is obtained from Theorem 14 with c = 0. 

This theorem helps explain the phenomenon of the efficiency of 
large trunk groups since even when a > x, B (x, a) is small as long as a 
is in a small neighborhood of x; thus, Theorem 14 shows that 
B(x, x + c Vz) ~ l/a yfx, x — * °o . 

Theorem 14 shows that the parameter c may be viewed as a standard- 
ized offered load measuring the deviation of a from x in units of V^- 
The value of this viewpoint is derived from the very simple approximat- 
ing form for B (x, a) ; thus, 

B(x, a)- 1 ~ aoV^ + ax + ~ (71) 

An application of this is to the computation of dB(x, a)/dx given in 
(92). Another advantage is the capability of computing B(x, a) by 
means of a single-entry table against the standardized offered load c 
rather than the usual double-entry table against x and a. 

Table I gives the values of a (c), ai(c), a 2 (c) for -3 ^ c ^ 4 in 
steps of 0.1 with the intention of covering a practical range of values 
of B(x, a). As an illustration, it is desired to compute fi(400, 378). 
Use of (71) with c = - 1.1 gives the result 0.0122 correct to the last 
figure. If c does not appear in the table, then interpolation is used. For 
example, to compute fi(400, 377.6) for which c = - 1.12 linear inter- 
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polation in the table of coefficients and use of (71) yields 0.0118 cor- 
rect to the last figure. The method, of course, is valid even when the 
number of trunks is nonintegral. Consider, for example, B (400.34, 
420) for which c = 0.98463. The result obtained by linear interpolation 
in the table is 0.0713 correct to half a unit of the last figure. 

The accuracy deteriorates when x is decreased or when c is large and 
negative. Thus, for 5(10, 8), one obtains 0.12144 as against the correct 
value 0.12166. In this case c = - 0.6325 is not too disadvantageous. 
The case 5(10, 5) for which c = - 1.58114 yields a much greater 
error, namely, 0.0256 as against the correct value 0.0184. This error 
occurs, however, for a small trunk group where exact calculation is 
quite feasible. To aid the delineation of suitable regions of (c, x) for 
which the table is accurate, a curve is given in Fig. 1 defining 5-percent 
error. When a computation is made from the table using any point 
(c, x) in the unbounded, unshaded region, the error incurred will be 
less than 5 percent of the true value of B(x, a). 

IV. DERIVATIVES AND INEQUALITIES 

It is desired to obtain formulas for the derivatives of B(z, a), with 
respect to z and a. 

Theorem 15: 

**(*>"> = I I - l + B(z, a)\ B(z, a), Rea > 0. 

da [a ) 

Proof: From (23), one has 

aB ^ a )~ 1 , f x e -««(i + u)'du -at* r-«(l + uYudu; (72) 
da Jo Jo 



hence, 



dB{z > a) ' = -B(z, a)- 1 - B(z + 1, a)- 1 + B(z, a)"'. (73) 
da a 



Use of Theorem 4 provides the relation 



dB(z, a)~ l z 



da a 



= - -B{z, a)- 1 -1+B(z, a)~K (74) 



Since 



agfe «)- m _ mz, a) (76) 

da da 



the result of the theorem follows from (74). 

ERLANG LOSS FUNCTION 541 



For the purpose of obtaining an approximate formula for the deriva- 
tive with respect to z, consider 

/(it) = aB(x, o)e-° u (l + u) x (76) 

in which a > 0, x > 0, and for which, by (23), 

J" f(u)du = 1. (77) 

It is convenient to introduce the random variable £ with density func- 
tion f(u). The power moments n r denned by 

H T = E?, r > (integral), (78) 

are given in the following theorem. 

Theorem 16: 

Mr = B(x, a) t (-l)'-< (Vj B(x + I, a)" 1 . 

Proof: Define a generating function 0(0 by 

0(0 = Ee't = aB(x, a) J* e-<°-' )u (l + u)*du, (79) 

then, since 

(a - t)B(x, a - t) f M e-<°-' )u (l + u)*du = 1, 



(80) 



one has 



aB(x, a) / Q i\ 

<K0 = 7 AP/ A" ^ 81) 

vw (a — 0B(s, a - t) 



Use of Theorem 10 in (81) provides the expansion 

*(i) = B(x, a)e- f *(? + T,ar* C (82) 

r=0 " ; 

Since 

♦«>-£§«', ( 83 ) 



(84) 



the coefficient of t r in the expansion of (82) in powers of t yields the 
required result. Thus 

Mr = *(*, a)r ! L y^Ty, f\ 

and the formula of the theorem follows. 
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Corollary : The central moments a r are given by 

a r = E(t - /xi) r = B(x, a) £ (-1)- 1 (]} fa + l)^B(x + I, a)~K 

Proof: The same as for Theorem 16 but considering the function 
e-«'0(O instead of <f>(t). 

An approximation to dB(x, a)/dx may now be obtained. 

Theorem 17: 

. , dB(x, a) , /t , v 1 a 2 

+ 5(if^-J M ' 0<9<1 ' 

Proof: From (23) and (76), one obtains 

-£(*, a)- 1 dB{ £ a) = E\n (1 + 5). (85) 

Since, by use of the mean value formula, 
In (1 + = In (1 + mi) + ^ (t - mi) - g 7jq^)i (t - Mx) 2 

+ H (i+ 1 Mi)« u " Ml)3 ~i gU ~ Ml)4 ' 0<d<1 ' (86) 

one has, from (85) and the corollary to Theorem 16, the required 
result. 

Corollary : 

D . . 35(a:, a) . . . 1 a 2 ,1 «3 

-b(«, «)- — jj- < in (i + n) - gjmu' + 3 ir+70 5 ' 

Proof: The error term of Theorem 17 is omitted. 

For ready reference the following formulas are given in which 
B = B(x, a), Bi = B(x + 1, a), B 2 = B(x + 2, a), B 3 = B(x + 3, a), 
£ 4 = B(x + 4, a). 

Mi = - 1 + SB," 1 , 

a 2 = (mi + l) 2 - 2(mi + lJBBf 1 + Bflr 1 , (87) 

a, = - (mi + l) 3 + 3 (mi + l)*BBr l - 3 (mi + l)BB2 l + BBf\ 
ai = (mi + l) 4 - 4fa + lYBBr 1 + 6(mi + l) 2 BB 2 - ! 

- 4 (mi + DBBf 1 + BBr 1 . 
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The evaluation of Br 1 , Bf 1 , B?\ Br 1 is facilitated by successive use 
of Theorem 4. 

An alternative method of obtaining dB(x, a)/dx is based on Theorem 

14. Let 

f(x, c) - B(x,a), a = x + cVx, (88) 

then, from Theorem 14, 
hence, 



g/fe C > ! ~ - f ^J aj { C )x-W\ X -> oo ; (89) 

3x ; =o 2 



ax 2Vx r w x 



(90) 



Thus, the computation of df(x, c)/dx is easily accomplished with the 
help of Table I and the formula 



a/^.- Mc) .^~-^U)-SM) 



(91) 



One now has 



dB(x, a) df(x, c) dB(x, a) / . ■ _g_\. ,q 2 \ 

~~ dx~ dx da V 2W 

A simple upper bound on -B(x, a)- l [dB(x, a)/dx~\ is given in the 
following theorem. 

Theorem 18 : 

_B(*,a)-*42fe^<in(l + ni). 

Proof: Since the function - In (1 + u) is convex for u ^ 0, the re- 
quired inequality follows from Jensen's inequality, namely, 

g(Ei) ^ Eg {!■) (93) 

valid for functions g(x) convex over the range of the random variable 
£, and (85). 

A function g(x) > is said to be log-convex over a set if In g(x) is 
convex over the set. It is known 19 that the sum of log-convex functions 
is log-convex and hence that the integral of a log-convex function with 
respect to a parameter is log-convex provided the function is log-convex 
for every value of the parameter. Since a necessary and sufficient 
condition that a twice-differentiable function be convex is the non- 
negativity of its second derivative over the corresponding set, one 
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derives the inequality 

g"g -g"^0 (94) 

as a necessary and sufficient condition that g > be log-convex. One 
now has 

Theorem 19: B(x, a) -1 , [_aB{x, a)~]~ l are log-convex junctions of x and 
of a, respectively, for a > and all x. 

Proof: The results are immediate from (23) and the observations that 
(1 + u) x is log-convex as a function of x for u ^ 0, and e~ au is log- 
convex as a function of a for u ^ 0. 

Corollary : 

d*B(x, a) \ dB(x,a) V 

n/ , [ n dB(x, a) . d 2 B{x, a) 1 . [ D , N , dB(x, a) l 2 
aB{Xi a) [2 — i^-i + a— £f-' J ^ [B(x, a) + a—£- J • 

Proo/: Use of (94). 

An immediate application of Theorem 19 is to the logarithmic inter- 
polation of B(x, a), that is, linear interpolation of In B(x, a). 

Theorem 20: Let a,b, p, q > 0, p + q — 1, then 

B(x, a)'B(y, a)" ^ B(px + qy, a), 
[aB(x, a)2'D>B(x, b)]« ^ (pa + go)B(x, pa + ob). 

Proo/: Jensen's inequality applied to —In B{x, a) and —In [aB(x, a)], 
respectively. 

An extension of the result of Theorem 20, for the purpose of obtain- 
ing an approximate formula for B (x, a) when x is not an integer, may 
be derived from the corollary to Theorem 16. Let N = [x], 8 = x — N, 
and a r be the central moments computed for the density function 

f(u) = aB(N, a)e-«"(l + u) N , (95) 

then one has 

Theorem 21 : 

B(x, a)" 1 = B(N, a)- 1 if ( * ) «r(l + m)'- + B(N, a)" 1 ( * ) aj, 

k even, |0| ^ 1. 

Proof: Let £ be the random variable with density function f(u), then 

B(x, a)- 1 = B(N, a)- l E(l + £) 4 . (96) 
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Since 

(1 + £) s = Z7 I ) (1 + Mi)*- p tt - Mi)' + ( I ) U - Mi)**, (97) 

the result follows from (96) and the corollary to Theorem 16. 
A useful special case of Theorem 21 is 

B (x, a) ^ ") 2 , (98) 

in which 
B = B(N,a), Bx = B(N+l,a), B 2 = B(N + 2,a). (99) 

V. CONCLUSION 

Further investigations would be desirable; for example, one would 
like to know the contour function g(z) for which B[z, g(z)~\ is constant. 
Truncation error formulas for the asymptotic expansions of Theorems 
13 and 14 would be useful ; also, the general structure of the coefficients 
a,j(c) of Theorem 14 should be determined. Asymptotic formulas of 
various types should be obtained for G,-(z, a) similar to those given for 
B{z, a). These formulas may then be used to obtain asymptotic results 
for its zeros which are needed in many transient and time-variable 
blocking analyses. 
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APPENDIX 

The function B(z, a) is related to the Poisson-Charlier poly- 
nomials 20-22 much used in telephone traffic studies. Let 

ft, (2, a) = 4t(z, a), 

*,(*,«) = a&*k«o, ( 10 °) 

then the Poisson-Charlier polynomials, Qj(z, a), are defined by 

fcfe a) = ffo «)(?/(«, a). (101) 
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The Taylor expansion 

*(*,« + «) = £-£*/(*«) (102) 

yields the generating function 23 

e-'(l+-V = E6A«)r,- (103) 

\ «/ y-o J * 

Thus, explicit formulas for G>(z, a) are 

■^^^(O'^O)' (104) 

The first few polynomials are 
G (z, a) = 1, 

Gi(z, a) = - [z - a), 

" (105) 

<?.(*, a) = i [z 2 - (2a + 1)2 + a 2 ], 

<?»(*, «) - - 3 L* 3 - 3(a + l)z 2 + (3a 2 + 3a + 2)z - a 3 ], 
a" 

A recurrence relation derived from (103) is 

G J+1 (z, a) = Z ~ J a ~ a Gj(z, «) - l Gj-i(z, a). (106) 

The polynomials, G>(z, a), possess many properties analogous to 
those of B(z, a)~ l . A Rodrigues formula is given in 



Theorem 22 : 



G j+k (z, a) = cT'e* ^s [e-WQfa a)]. 



Proof: One has from (100) 



and hence 



*,+*(*,«) -g^ */(•,«), (107) 



<W*, «)*& a) = ^ [#(«, «)<?,(*, a)]. (108) 



The result follows on use of (24). 
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Corollary : 

G,(z, a) = cT'ef — [r^Wj 

Proo/: 

G («, a) = 1. (109) 

An integral representation for (?y(z, a) is given in 
Theorem 28 : 

g } (-z, a ) = (-*)' jSg J* «- av (i + yyy-'dy, 

Re a > 0, Re 2 > 0. 
Proof: From (20), one has 

rv = =^r [" er'^y-idy. (110) 

l^J ./o 

Substitution of (110) into the corollary of Theorem 22 yields the result. 
Theorem 23 permits obtaining a Wittaker function representation. 

Theorem 24 : 

Gi(z,a) = (-l)4rW H "W tfHt i,/i I (H)/i(«). 

Proo/: Comparison of Theorem 23 with (32) and replacement of —2 
by 2. 

The representation of Theorem 24 remains valid, by analytic con- 
tinuation, even outside the region of convergence of the integral of (32). 

Corollary 1 : 

B(N,a)~ l - (-l) N G N (-l, a). 

Proof: Comparison of Theorems 7 and 24. 
Corollary 2 : 

Gj{Z, a) = ^Gj-xiz - 1, a) - (?y_i(2, a). 

Proo/: Substitution of the representation of Theorem 24 into the re- 
currence relation 24 

W k , m (z) = ^W k - hm -i(z) + (*-*; + m)W k ^, m {z) (111) 
yields the result. 
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Corollary 3 : 

Gj(z, a) = Gj(z - 1, «) + |(?,-i(* - l, a ). 

Proof: Same as for Corollary 2, except the following recurrence rela- 
tion is used : 

W k .M = 4~zW k - h . m+h {z) + (| - k - m)Wu- X , m {z). (112) 
A representation of G,{z, a) in terms of B(z, a) is given in 

Theorem 25: 

<?„(-;,«) = (-^"u^^^h)- 

Proof: From (23), 

__L = T e -a y(1 + y)N d y (113) 

aB{N, a) Jo 
and Theorem 23, 

<?#(->, «) = (-1)* ( - "\ } , / o " r-»(l + yVy^dy, (114) 

one has the result on use of 

fjj-l roa fee 

|— j j o e-«v(l + y )'^ = (-l)'-i I e-""(l + y) N y^dy. (115) 

The Poisson-Charlier polynomials possess addition formulas similar 
to those of B(z, a) -1 as given in Theorems 9 and 10. 

Theorem 26: 

e,k« + = (i + t)- J 'f o (^ G/ .,(,, «)(-i).(i)'. 

Proof: Use of (103) shows that the system of functions 

has the generating function e at (l — t)*, and hence 25 form an Appell 
system with respect to a, thus, 



i [ Sr "*'<* a) ] = b^W\ ai ~ lGi ~ l{z ' a) - (1 



16) 



The Taylor expansion of [(— l)1/j I] (a + t)'Gj{z, a, + t) in powers of 
/ now yields the required result. 
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(?,.(*, « + t) = (l + *- )"V £ <?*,(*, «) f. 



77ie series is permanently convergent. 

Proof: Use of (107) and Taylor's expansion yields 

v =o v \ 
The result is now obtained from (101) and (117). Since, from (104), 

OM4~(z)', i-«. ( 118 > 

the convergence is permanent. 

An asymptotic expansion is given by 



Theorem 28 : 

Gj(z,a)~ i-iy 



1+ t (-D 



jMjM 



=0 " i« 



l«» 



a — » oo. 



Proof: The result is obtained directly from (104). It may also be ob- 
tained from Theorem 24 and (45). 
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